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Ion traps are a versatile tool to study nonequilibrium statistical physics, due to the tunability of 
dissipation and nonlinearity. We propose an experiment with a chain of trapped ions, where dissipa- 
tion is provided by laser heating and cooling, while nonlinearity is provided by trap anharmonicity 
and beam shaping. The collective dynamics are governed by an equation similar to the complex 
Ginzburg-Landau equation, except that the reactive nature of the coupling leads to qualitatively 
different behavior. The system has the unusual feature of being both oscillatory and excitable at 
the same time. We account for noise from spontaneous emission and find that the patterns are 
observable for realistic experimental parameters. Our scheme also allows controllable experiments 
with noise and quenched disorder. 

PACS numbers: 



Pattern formation is the emergence of structure in a 
nonlinear medium far from equilibrium [J, Q . This phe- 
nomenon occurs in many settings, including fluids, chem- 
ical reactions, plasmas, and biological tissues. In tra- 
ditional pattern-forming systems, the collective behav- 
ior is set by the material properties, and theoretical de- 
scriptions are often phenomenological. For example, in 
the Belousov-Zhabotinsky reaction, the concentrations of 
chemical reagents oscillate in time and produce traveling 
waves. It is a complicated reaction involving many inter- 
mediate states and rate constants. Hence, it is difficult 
to experimentally control the behavior of the system. 

On the other hand, ion traps allow an unprecedented 
level of control using optical and electrostatic forces and 
have led to impressive experiments in quantum comput- 
ing Q and quantum simulation In this paper, we 
show how ions are also useful for studying pattern for- 
mation. The advantage of using ions here is the ability to 
tune dissipation and nonlinearity in situ, thereby having 
more experimental control and being able to see different 
effects within the same system. 

We show how patterns arise in a chain of ions driven far 
from equilibrium. The collective dynamics are governed 
by an equation similar to the complex Ginzburg-Landau 
equation, which is one of the most studied equations in 
physics. However, the presence of only reactive coupling 
in the ion chain leads to novel behavior: while other sys- 
tems are either oscillatory or excitable, the ion chain can 
be both at the same time. Our scheme is also useful 
for studying synchronization and Anderson localization. 
Our work is motivated by recent experiments on the non- 
linear dynamics of single ions: a phonon laser and a 
Duffing oscillator [7|. We note that the model described 
in this paper is also applicable to an array of nanome- 
chanical resonators 

First we describe the proposed experimental setup. A 
linear Paul trap uses an RF electric field for radial con- 
finement and a DC field for axial confinement We 
use a segmented trap, which has many DC electrodes in 




FIG. 1: Chain of ions, each damped by a red-detuned beam 
and excited by a blue-detuned beam. The blue beam is along 
the ion chain, while the red beams are at an angle. The 
intensity of the red beams is lowest at the trap center. Each 
ion is in its own anharmonic DC potential well (not shown). 



order to create many trapping regions [lOj], and thereby 
make a chain of ions, each in its own potential well. By 
changing the DC voltages, one can tune the shape of the 
potential for each ion and thus control the nonlinearity 
[III [l2| . Let x be the axial displacement of an ion from 
its trap center, d the distance between trap centers, lo 
the harmonic trap frequency, and a the coefficient of the 
anharmonic quartic term in the trap potential. 

We apply near-resonant laser beams to heat and cool 
each ion (Fig. [T}. When the laser frequency is above 
resonance (blue-detuned), the ion feels an anti-damping 
force due to the Doppler effect. When the laser frequency 
is below resonance (red-detuned), the ion feels a damping 
force [l3(. We apply a blue-detuned beam in the axial 
direction along the ion chain. For each ion, we also apply 
a red-detuned beam at an angle <f> with respect to the 
trap axis. The red beam is shaped so that the intensity 
is lowest at the trap center. The ion is then heated near 
the center and increasingly cooled away from the center, 
so the ion oscillates with a large amplitude, determined 
by the balance of heating and cooling. The dissipation is 
easily tuned by changing the beam intensity. 

We use a singly- ionized, two- level atom of mass m with 
a dipole transition of wavelength A = 2n/k and linewidth 
7. Each red and blue beam has detuning —Auj and +Aw 
and intensity Jr and Ib, respectively. Let I s be the sat- 



2 



uration intensity. All beams have the same polarization. 

We use low laser intensities so that the ion is unsatu- 
rated. We assume the radial motion is cooled near the 
Doppler limit (due to the projection of the red beams), 
so that the Doppler shift is due mainly to the axial mo- 
tion. For mathematical convenience, we use counter- 
propagating beams (but in practice one would use sin- 
gle beams). Then the total optical force on an ion is 
calculated rigorously to be [H, 



( 7 2 + 4Aw 2 ) 2 



Ir 2 
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under the additional assumptions k\x\ -C 7 2 /4 and 
k\x\ <C A.lu. To get position-dependent damping, we 
choose to vary Ir quadratically along the trap axis, 



£cos<fi 
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where I is the characteristic length of the intensity gradi- 
ent. The intensity profile does not need to take this form 
or even be symmetric. In fact, a different profile may 
lead to interesting higher-order terms in Eq. ([5]) |14j . 

The ions are coupled through Coulomb repulsion. If 
the displacements are small relative to the inter-ion dis- 
tance [£ <C d), the interaction is linear, and the free-space 
coupling decreases with the cube of the distance. Numer- 
ically, we find that interactions farther than the nearest 
neighbor do not affect the overall dynamics much, so we 
assume only nearest-neighbor interactions. The equa- 
tions of motion are, 
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n = 1, . . . , N, where k e is the Coulomb constant, e is the 
proton charge, fj, is the damping coefficient, 
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and Xn(t) is the noise. In this scheme, the inherent 
source of noise is spontaneous emission, since each emis- 
sion causes a momentum kick hk in a random direction. 

We work in the regime where the nonlinearities and 
interactions are small perturbations to the harmonic mo- 
tion. We write x n (t) = 2£ ReL4 n (t)e _i ""], so that the 
complex amplitude A n encodes the slowly varying am- 
plitude and phase of the underlying harmonic oscilla- 
tions. In the Supplemental Material, we find the am- 
plitude equation, 



dA„ 
dt 



A n + ib{-2A n + An-! + A n+1 ) - (1 + ic)\A n \ 2 A r 



-fo«(t) + ia«{t)]A n + kf (f) + «r*(i)] 
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where t — /it/2 is rescaled time, b is the coupling, and c 
relates how an ion's amplitude affects its frequency. We 
stress that b and c are directly related to experimental 
settings. In the absence of coupling, each ion oscillates 
with amplitude \A\ = 1, which corresponds to an ampli- 
tude of 2£ in x. The noise functions are due to sponta- 
neous emission and represent scattering by the red beams 
(r] R ,a R ) and blue beams (rj B ,<j b ), 
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where H is a dimensionless measure of the noise. 

We now examine the spatiotemporal properties of 
Eq. ([5]), ignoring the effect of noise for now. First note 
that the equation is symmetric under the transformation 
b,c,A n —} —b,—c,A n . In the continuum limit, we let 
A n -> A(X) and 



dA 
~~dl 



A + ib 



d 2 A 
dX 2 
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This is similar to the complex Ginzburg-Landau equation 
(CGLE) except that the coefficient of d 2 A/dX 2 is 
purely imaginary. This is because the Coulomb force is 
reactive, while the CGLE includes both reactive and dis- 
sipative interactions. This greatly modifies the behavior, 
as seen below. 

According to Eq. (fTUj) . a plane wave solution A(X, t) = 
Fe i(QX-wi) sat i snes F =1 and uj = bQ 2 + c. By lineariz- 
ing around this solution Q, we find that the condition 
for stability is be > 0, and as long as this is fulfilled, there 
is no restriction on the wave number Q. Since Eq. (jf 01) is 
in the continuum limit, we expect long wavelength waves 
(wavelength at least several ions) when be > 0. 

When be < 0, it turns out that Eq. (O allows short 
wavelength waves that are not captured in the continuum 
limit. Define A n = —A n for even n and A n = A* for odd 
n and consider the continuum limit of the transformed 
system. A plane wave A(X) = Fe 1 ^*'^ satisfies F = 
1 and Cj = b(Q 2 — 4) — c, and the stability condition for 
long- wavelength waves is be < 0. A long wavelength wave 
in A corresponds to a very short wavelength wave in A. 

Therefore, Eq. ([5|) has stable plane waves for all values 
of b and c: long wavelength for be > and short wave- 
length for be < 0. This behavior is different from the 
CGLE, which has stable plane waves for be > — 1 and is 
otherwise chaotic 15|- Another difference is that in the 



CGLE, only a band of Q is stable. 

However, boundary conditions affect the selection of 
plane waves. With periodic boundary conditions, the 



3 




20 30 
oscillator index 



ial 



:-> 




390 _ 400 

rescaled time t 



FIG. 3: System of two oscillators with 6 = 1. (a) Bifurcation 
diagram as c varies. There are supercritical pitchfork bifurca- 
tions at (c, A(9) = ( — 1,0) and (1, 7r). There are supercritical 
Hopf bifurcations at (1.22,2.82), (1.22,3.47), (-1.22,-0.32), 
and ( — 1.22,0.32), which give rise to stable limit cycles (not 
shown). Solid and dashed lines denote stable and unstable 
fixed points, respectively, (b) A limit cycle at c = 1.4. 



FIG. 2: Spacetime plot of chain of 50 ions with nearest neigh- 
bor interactions and open boundaries. Re A is plotted using 
the color scale in the side bar. A is the complex amplitude 
of the underlying harmonic oscillations, (a) b = 1 and c = 1 
showing uniform phase synchrony, (b) Same, but with the 
expected noise from spontaneous emission, (c) b — 1 and 
c = — 1 showing anti-phase structure. 



wave number (Q or Q) is a multiple of 2-k/N. Open 
boundary conditions are simpler to implement and are 
equivalent to setting 44 = at the boundaries. Thus 
when be > 0, the only allowed plane wave is the Q = 
wave, in which the ions are uniformly in-phase (Fig. [2^,). 
(This also occurs in the Belousov-Zhabotinsky reaction 
in the absence of spatial inhomogeneities Q). One can 
induce Q ^ waves by, for example, changing the trap 
frequency uj q of one ion. The be < case is different, 
because open boundary conditions mean setting A = at 
the boundaries. The final state is not a pure plane wave 
but a more complicated structure, in which the ions are 
almost uniformly anti-phase (Fig. &p). 

We now examine the dynamics of two coupled ions, 
expanding on previous work [3, 17 1. First write Eq. ([5]) 



in terms of the real amplitudes r\ , T2 and phase difference 
A9 = 9 2 - 6»i, where A n = r n e~ ie " , 
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This system is symmetric under the following transforma- 
tions: { T-i, r-2, AO -> r 2 ,ri,-A0}, {c, A9 -> -c,tt- A9}, 
and {b, A9 — > —b,n + A9}. There are fixed points at 
(AS, 7-1,7-2) = (0,1,1) and (tt, 1, 1), corresponding to in- 
phase and anti-phase motion. There is another set of 
fixed points that correspond to roots of a quartic poly- 
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i-(c 2 -l)(c 2 -3K-( C 2 
(6 2 + l)(c 2 + l)r 2 + 6 2 (6 2 + l) , 
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which may be solved numerically. The bifurcation di- 
agram is quite rich: as b and c change, saddle-node, 
pitchfork, and Hopf bifurcations appear, disappear, and 
change criticality. An example is shown in Fig. For 
some values of b and c, there are supercritical Hopf bi- 
furcations to stable limit cycles, in which the amplitudes 
and relative phase oscillate (Fig. [3b) . The system is at 
least bistable for \c\ < \b\, but certain values of b and c 
have four stable fixed points. 

As N increases, there are still in-phase (A9 — 0) and 
anti-phase (A9 w tt) fixed points, although the region of 
multistability in be space shrinks. For |c| < |6|, there are 
also multiple stable limit cycles, in which the entire chain 
has the same average frequency or the chain is divided 
into regions of different frequencies. 

A large ion chain is excitable in a novel way. An ex- 
citable medium has the property that the uniform state 
is stable to weak perturbations, but a perturbation that 
exceeds a threshold grows rapidly and then decays. Usu- 
ally, a medium is either oscillatory or excitable 2] . How- 
ever, the ion chain is both oscillatory and excitable at 
the same time. Suppose that be > and the chain is in 
the Q = state (in-phase), with large. We perturb 
the end ion A\ by SA. The Q — state is linearly sta- 
ble so a small perturbation decays away. But if 8 A is 
greater than a threshold, it generates a localized pulse of 
anti-phase oscillations (Fig. 01 . The pulse travels across 
the system, bouncing off the boundaries until it decays. 
This type of excitability differs from traditional examples 
(like neurons and heart tissue) , because the pulse is made 
of an alternating phase structure instead of an increased 
chemical concentration. 

The excitability can be intuitively understood from the 
fact that when be > and |c| < both in-phase and 
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FIG. 4: Spacetime plot of Re A for chain of 50 ions showing 
an excitation pulse for b — 1 and c = 0.2. Perturbing the first 
oscillator beyond a threshold generates a pulse of anti-phase 
oscillation. The initial conditions were Ai = — 1 and A n = 1 
for the rest. The right panel is a zoomed-in view. Color scale 
is the same as in Fig. [2] 



fluctuate not only in time but in space as well. In our 
scheme, the noise for each ion is independent and may 
be easily tuned by changing the detuning Auj in Eq. © . 
One could observe, for example, noise-induced transitions 
between stable fixed points in a small chain. 

Another interesting use of the tunability is to study 
the effect of quenched disorder. A previous work studied 
the mean-field version of Eq. ([5]) with random harmonic 
frequencies uj and found that as b and c change, the sys- 
tem undergoes continuous and discontinuous phase tran- 
sitions between the unsynchronized and synchronized 
states [H, . It would be interesting to study the lower- 
dimensional versions. We note that synchronization of 
disparate oscillators is an important topic throughout 
science 2^, 21 1. Furthermore, when the variance of lo is 
small, Eq. (|5|) may be mapped via the Cole-Hopf trans- 
formation to the Schrddinger equation for a particle in a 
random potential [22|, HH • The resulting pattern forma- 
tion reflects the phenomenon of Anderson localization. 

We thank H. Haffner for helpful comments. This work 
was supported by Boeing and NSF grant DMR-1003337. 



anti-phase oscillations may be stable for small chains, 
while only in-phase oscillation is stable for large chains. 
Thus, a local region within a large chain may be anti- 
phase for a short amount of time. A mathematical de- 
scription of this phenomenon is left open for future work. 

An important question is whether the patterns de- 
scribed above would be observable for realistic experi- 
mental settings that satisfy all the theoretical assump- 
tions. Indeed, we find that the patterns are visible above 
the noise from spontaneous emission. For example, the 
ion 24 Mg + has an SV/2 — P3/2 dipole transition of wave- 
length A = 279.6 nm and linewidth 7/2tt = 42 MHz. 
Letting uj /2tt = 100 kHz, I B /h = 0.05, Aw = 67, 
4> = tt/4, I = 30 ^m, a /47r 2 = 10 15 Hz 2 /m 2 , and 
d = 500 /im, one finds b = 1.0, c = 1.1, and H = 5x 10~ 4 . 
Figure [2Jd shows that the Q = state is clearly vis- 
ible above the noise. Also, since many experiments 
are already using large-scale traps for quantum informa- 
tion (lpj . it should be straightforward to implement our 
scheme with many ions. 

In the experiment, one can measure the amplitude and 
phase of A n by recording when and where the ions scatter 
photons [H, 0|- Changes in A n occur on a time scale 
t ~ 1/m> which is much slower than that of the harmonic 
oscillation (~ 1/uj ). Thus, one can observe dynamical 
effects, such as limit cycles and excitation pulses. The 
entire be space can be explored by, for example, tuning 
the parameters Ib and a Q . 

It would be interesting to turn up the noise to see what 
happens. It is known that adding noise to a sp atially ex- 
tended system may have nontrivial effects [18J. However, 
it is usually experimentally difficult to make the noise 
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Supplemental Material 

This appendix provides supporting calculations for the main text. In the first section, we derive the amplitude 
equation using perturbation theory. In the second section, we calculate the noise from spontaneous emission. 



Derivation of amplitude equation 



Here we use perturbation theory to calculate how the weak nonlinearities and interactions in Eq. (3) of the main 
text change the amplitude and phase of the harmonic oscillations on a long time scale. We use the method of averaging 
with amplitude and phase variables because the noise from spontaneous emission depends on the amplitude. 

First we rescale Eq. (3) of the main text with r = uj t and y n = x n /£, 



d 2 d 
= -r^Vn + Vn + ay n - v(l - y n )—y n + D[(y n - y n _ x ) + (y n - y n+ i)} + C„(t) 



(15) 



where a — a £ 2 /ui 2 , v = fj,/cJ > D — 2k e e 2 /md 3 uj 2 , and Cn( T ) is the noise. Equation ([15)) describes a chain of van der 
Pol-Duffing oscillators. 

Let y„ — r' n (r) cos [r + # n (r)], where r' n and 9 n change slowly (r' « r', f' « f', 9 < 1 J « fl). Substituting this 
into Eq. (|15l) and keeping the leading terms, 

= -2f' n sin(r + 6 n ) - 2r'J n cos(t + 9 n ) + ar' 3 cos 3 (t + 6 n ) + v[l - r' 2 cos 2 (r + 8 n )]r' n sin(r + 9 n ) 

+D[2r' n cos(t + 9 n ) - r' n _ 1 cos(r + 6 n -x) ~ r' n+1 cos(r + 6 n+ i)} + C„(r) . (16) 

Then multiply each equation by sin(r + n ) and integrate over the time interval [r, r + At], where At is a multiple 
of 2tt, 



dr 



,.'2 



mdi((? w _i - 6 n ) + r' n+1 sm(6 n+1 - 6 n )] + ^(t) 
Then multiply each equation by cos(t + 9 n ) and integrate similarly, 



dK 
dr 



3 ,2 D 
-ar H 

8 " 2 
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These equations describe how r' n and 9 n evolve. The noise functions are, 

0-> = -T- / r+AT rfT'C n (T')sin(T' + ^) 



1 
A7 
1 



T+Al 



dT , Cn(T , )cOs(T / +9 n ) 



r' n A r 

To put Eqs. (|17l) and (fl8l) in simpler form, rescale time (t = vt/2) and amplitude (r n — r' n /2), 
dt 

d9 n 2 



dt 



(1 - r n )r n + b[r n -i sin(6» n _i 
b 



+ r n+ i sin(^ n+ i 



cr 



2 cos( 



\ r n+l , 
n-l) C0S( 



)]+V n {t) 



(17) 

(18) 

(19) 
(20) 

(21) 
(22) 



where b = D/v, c = 3a/ is, and ip^ and tpf^ are the rescaled noise functions. Then write everything in terms of a 
complex amplitude A n = r n e~ l9n , so that y n {T~) = 2Re[A n (r)e~ lT }. (We use e~ t6 " instead of e l6n in order to match 
up with the sign convention in Ref. (la|-) Then A n (t) evolves according to, 

dA n 



dt 



A n + ib{~2A n + A n _i + A n+1 ) - (1 + ic)\A n \ 2 A n + ip£(t, A n ) , 



(23) 



where ip^ is the complex-valued noise function. 
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Noise from spontaneous emission 



Here we calculate the expected noise from spontaneous emission. When an ion absorbs a photon from a laser, it 
gets a momentum kick in the direction of the laser, and when it spontaneously emits the photon, it gets a momentum 
kick in a random direction. Spontaneous emission is the inherent source of noise in our scheme, so we explain how to 
represent it with the noise term ip A in Eq. ([23]). 

There are two factors that must be taken into account. First, for the experimental conditions assumed in the text, 
an ion scatters on the order of one photon per oscillation cycle. Thus, the noise is a sequence of occasional impulses 
happening at random times. Second, the noise is position dependent due to the intensity gradient of the red beams. 

We just consider a single ion, since the noise for each ion is independent and identically distributed. Each scattering 
event happens at a random time, and the spontaneous emission of a photon causes a momentum kick hk in a random 
direction. Suppose the ion scatters photons at times t n . Then the noise in Eq. (3) of the main text is 

X(t) = -$>(*-*n)9n, (24) 
n 

where q n is a random variable (with variance a q ) for the projection of a momentum kick along the trap axis. Each 
kick is independent {{qjqk) = $jk)- F° r simplicity, we assume that the emission is isotropic (cr 2 = 1/3), although there 
is a slight anisotropy relative to the laser direction [24| . 

With the assumptions on experimental parameters given in the text, the scattering rate T may be calculated 
rigorously from the Optical Bloch Equations [U [25| , 



T(x) 



7 3 f Ir(x) I b 



7 2 + 4Awf, 7 2 



(25) 



The first and second terms correspond to scattering by the red and blue beams, respectively. Note that T depends 
on position and is independent of velocity to first order. 

After rescaling (r = uj Q t and y = x/i) to get Eq. (fT5"]l . the noise is 



hk x -> 



and the scattering rate becomes, 



CM = ^j^-r n)qn , (26) 



r(y) = r R ( y ) + r B (27) 

fs = if 72+^ ' (29) 



where we have used the intensity relation given in Eq. (2) of the main text. 
To calculate the amplitude noise £ r , plug Eq. (|2"B"j) into Eq. (pi 



= a J E <7nSin(T„ + 0), (30) 

t<t„<t+At 

where r„ is the time of a scattering event. Since the damping is weak, the ion scatters on the order of one photon in 
an oscillation cycle (r ~ lu /2it), so there is significant time between scattering events. This means that the phase of 
oscillation at which a scattering event occurs is approximately uncorrelated with the phase of the next event. Each 
scattering event has a random projection and phase. Thus, the sum in Eq. (|30[) is over independent samples of the 
random variable w n = q n u n , where u n = sinr n (ignoring the unimportant phase offset 9 for now). 

Now we find io n 's distribution p w . The intensity gradient of the red beams causes them to scatter more at certain 
phases within a cycle, while the blue beams scatter uniformly. Thus, p w is actually a weighted average of red and 
blue components. First we find the distribution of scattering times t„ (mod 2-7r) from the intensity profiles, 

i cos 2 r red 

cat) ; , (31) 

77- blue 

At: 
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since y = r' cost. Thus the distribution of u n = sinr„ is 

Pu{u) = p 



dr 
du 



7T 



red 



blue 



(32) 
(33) 



for \u\ < 1. Since we assume isotropic spontaneous emission, the distribution of the projection q n is p q (q) = 1/2 for 
\q\ < 1. Then the distribution of w n = q n u n is 



Pw(w) = du dq p u (u)pg(q)6(w - uq) 



Ilogi+^IH 

7T & | W | 



red 
blue 



(34) 
(35) 



for \w\ < 1. The variance of w„ is 



12 



red 
blue 



(36) 



To find the phase noise £ e , plug Eq. (|26]) into Eq. ([20]), 



^2 In cos(r„ + 6») 



(37) 



r<r I1 <r+Ar 

and go through the same process to find the variance of v n = q n cosr„, 

red 



i blue 

6 



(38) 



Although w n and v n come from the same scattering event, they are statistically uncorrelatcd because (sinr n cosr„) = 
0. 

We let the time interval Ar be large enough to include many scattering events but smaller than the characteristic 
time scales in Eqs. (fTTl) and (fT8|). We average the scattering rate T(y) in Eq. ([27)) over At to find the time-averaged 
scattering rates of the red beams (Tr) and blue beams (Tb), 



1 



Ib 
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T R = 



2oj \coscj)J I s 7 2 + 4 Aw 2 
1 Ib 7 3 



uj I s 7 2 + 4Aw 2 



(39) 
(40) 



Tr depends on r' due to the intensity gradient of the red beam. Then the amplitude and phase noises are Gaussian 
and described by, 



(fWf(r')) = 

(£V)£V)> - 



hk 
moj £ 
1 / hk 
r' 2 V maiJ 



1- 1- 

4 r * + 6^ 



5{t-t') 



(41) 
(42) 



They are uncorrelated with each other: (£, r (t)£, 6 (t 1 )) = 0. 
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After rescaling (t — vt/2, r — r'/2) to get Eqs. (j2Tj) and ([221) . the noises become, 



I f hk \ 2 f I - I 



<V r ( W)> = 5- 7 h^f fl + -f B *(t - F) (43) 

2^ V mu) t J \ 12 6 



1 f hk \ 2 1, 



<VW(F)) = 5-2 7 / R + -f B «(* - F) . (44) 

2^7^ \mu) t J \4 6 / 

Again, (ip r (F)^ 6 (F)) — 0. Finally, the complex-valued noise in Eq. (1231) is, 

V A M) = b ? K © + z ( T fl (t)]A+[r ? B (<)+ i a B (t)] (45) 

fo R (t> fl (F)) = -4-7<K? -F) (46) 

COS^ 

(* R (t)* R (i')) = -^r-S(t-F) (47) 
cos -2 <p 

(v B (t) V B (i')) = (<J B (F)a B (F))=H6(t-F) (48) 
ft(7 2 +4A^) 

where H is a measure of the noise, and we have simplified using Eq. (4) of the main text. The noise functions for 
the red beams (rj R ,a R ) and blue beams (j] B ,a B ) are all uncorrelated with each other. The noise from the red beams 
increases with amplitude and causes more phase noise than amplitude noise. 



